Evidence of the virtual Anderson transition in a narrow impurity band of 
p-GaAs/AlGaAs quantum wells: £4 conductivity and electric breakdown at low 

temperatures 
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In highly doped uncompensated p-type layers within the central part of GaAs/AlGaAs quantum 
wells at low temperatures we observed an activated behavior of the conductivity with low activation 
energies (1-3) meV which can not be ascribed to standard mechanisms. We attribute this behavior 
to the delocalization of hole states near the maximum of the narrow impurity band in the sense of 
the Anderson transition. Low temperature conduction £4 is supported by an activation of minority 
carriers - electrons (resulting from a weak compensation by back-ground defects) - from the Fermi 
level to the band of delocalized states mentioned above. The corresponding behavior can be specified 
as virtual Anderson transition. Low temperature transport (< 4 K) exhibits also strong nonlinearity 
of a breakdown type characterized in particular by S-shaped I-V curve. The nonlinearity is observed 
in unexpectedly low fields (< 10 V/cm). Such a behavior can be explained by a simple model 
implying an impact ionization of the localized states of the minority carriers mentioned above to 
the band of Anderson-delocalized states. 

PACS numbers: 73.21.-b 



I. INTRODUCTION 

The problem of a conductivity via localized states 
was addressed during a long time; however it mostly 
concerned a hopping conductivity. At the same time 
experimen1>i)^)^ demonstrated that a transition to metal- 
lic conductivity (or, in the case of 2D structures, to 
weakly localized conductivity) can occur within the im- 
purity band still split from the conductance or valence 
band. Recent observation of superconductivity in B- 
doped diamond which was supported by the states of 
impurity band^ once more stimulated an interest to the 
conductivity over localized states within strongly doped 
semiconductors . 

Note that typically the situation of the compensated 
samples was studied where the Fermi level was situated 
in a region of large density of states. Metal-insulator 
transition in this case is associated with Anderson tran- 
sition or, at least, with Mott-Anderson transition. As 
for the weakly compensated materials, until now a situ- 
ation of the metal-insulator transition is not completely 
clear. Indeed, without a compensation a finite conductiv- 
ity (including a hopping regime) can exist only via the 
states of the upper Hubbard band and in this case we 
deal with a "pure" Mott transition. However a presence 
of even small number of compensating defects allows a 
finite hopping conductivity over single-occupied states. 
Indeed, in this case a few "holes" within the impurity 
band, created by the compensating defects, are activated 
from the Fermi level to the maximum of the impurity 
band (nearest neighbor hopping). At the same time one 
can expect that if the concentration of dopants is large 
enough, then the carrier states within the impurity band 
become to be delocalized. Although the transport of the 
"majority carriers" is still forbidden due to strong on- 
site correlations, the activated "minority carriers" (aris- 



ing from the compensating defects) can propagate as it 
would be in a metal. So we deal with "virtual" Anderson 
transition: it happens in a situation where the delocalized 
states are still above the chemical potential, but the de- 
localized states do not exist until the concentration does 
not exceed some critical one. Since the disorder potential 
within the weakly compensated samples is weak enough, 
it allows to expect that such a "virtual" Anderson tran- 
sition can take place for the dopant concentrations suf- 
ficiently lower than for the Anderson transition for the 
"majority carriers". It is important that for 2D struc- 
tures the compensating defects can be situated outside 
of the 2D layer which also decrease the disorder poten- 
tial. 

In the case of n-doped materials the "majority carriers" 
within the impurity band are associated with electrons 
while the "minority carriers" are associated with empty 
donors ("holes"). At the same time for p-doped mate- 
rials the "majority carriers" correspond to holes while 
the "minority carriers" are related to negative acceptors, 
that is to electrons. It is important that the effective 
mass of the minority carriers is expected to be positive. 
Indeed, for the "majority carriers" the impurity band is 
nearly completely occupied while the " minority carriers" 
are situated close to the top of the corresponding energy 
band where the effective mass of the majority carriers is 
negative. Correspondingly, for the minority carriers the 
effective mass has the opposite sign with respect to the 
one of the majority carriers. 

The scenario of the metal-insulator transition (of Mott 
type or of Anderson type) depends in real systems on 
the specific parameters of the material (compensation de- 
gree and related scale of the disorder potential, effective 
masses of the carriers etc. In general, the problem now 
is far from a complete understanding which makes corre- 
sponding studies to be actual. 

In our recent publication- we reported an observa- 
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tion of activated conductivity at low temperatures in Be- 
doped uncompensated GaAs/AlGaAs single and multiple 
quantum wells with unusually small activation energies 
(more than an order of magnitude smaller than the Bohr 
energy of the dopant). It was shown that this tempera- 
ture behavior can not be associated with nearest neigh- 
bor hopping (£3 conductivity). We attributed such a be- 
havior to the "virtual" Anderson transition mentioned 
above. Namely, we believe that the observed activation 
energy corresponds to activation of the "minority carri- 
ers" from the Fermi level to delocalized states near the 
maximum of the impurity band. In the present paper 
we give detailed arguments in support of our conclusion 
including our new experimental results on low tempera- 
ture low field breakdown effects observed in the samples 
demonstrating the behavior mentioned above. To the 
best of our knowledge, the breakdown behavior and S- 
shaped 1-V curves observed in our experiments can not 
exist in the hopping regime. At the same time the scale 
of the electric fields (less than 10 V/cm) does not allow 
to relate these effects to breakdown to the valence band 
which would need much stronger electric fields. 



II. EXPERIMENT 

The technique of growth of multilayered structures 
with a help of MBE method was described in our paper—. 
The structures contained 1,5,20 GaAs quantum wells 
with widths 15 nm separated by barriers of Alo.3Gao.7As 
with widths 100 nm. The middle region of the wells was 
doped by p-type impurities (Be), the volume impurity 
concentration was controlled during the growth and var- 
ied from 1 • IOI8 atoms/cm"^ up to 2 • lO^^ atoms/cm'^ 
(Table 1). The critical concentration for a bulk p-type 
GaAS is 2 • lO^^ cm~^, i.e. the concentrations mentioned 
above are of the order or some less than the critical one. 
The compensation degree K = Nd/Na < 0,01 was suf- 
ficiently small and supposedly controlled by defects sit- 
uated at the edges of the quantum wells and within the 
barriers. All the samples parameters are given in Table 
1. The column Na gives the bulk acceptor concentration 
which was controlled during epitaxial doping. 



the conductivity and Hall effect with an activation en- 
ergy El ~ 16 — 20 meV (intermediate between Eq and 
eo/2, where eq ~ 28meV - is an energy of isolated accep- 
tor in GaAs). Note that in the samples suffering Mott 
transition, e^? is situated between the centers of the up- 
per and lower Hubbard bands: ef — £q — U/2, where U 
- is the Hubbard energy which for 2D is lower than Eq. 
Earlier we estimated the binding energy of doubly occu- 
pied state for the wells with a width 15 nm as 10 meV— , 
i.e. U = 18 meV; it gives Ei? ~ 19 meV, which agrees 
with observed values of £1 . 

For the rest of the samples with concentration of ac- 
ceptors close but still less than a critical one we observed 
a pronounced activated behavior of the conductivity at 
temperatures (10-1, 3K) characterized by a low activation 
energy £4 ~ 1 — 3 meV (Fig.l). 
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FIG. 1: Temperature dependences of resistivity for samples 
with a constant activation energy and for sample no. 200N2 
that is in the weak localization mode [3]. 

It was most pronounced for the single well structures. 
The value of £4 decreased with an increase of acceptor 
concentration (samples 945,946 ). As it is seen from Fig. 
1 and from the Table, the value of the preexponential 
CTo was only by a factor 10 - 30 lower than the universal 
quantum limit for 2 D : ctq < e^/Zi = 4 ■ 10~^O/im~i. 
Fig. 2 demonstrates temperature behavior of the Hall 
mobility for these samples. 



N 


number of wells 


well width, nm 


P300K,cm 2 


NA,cm ^ 


£i,meV 


£4,meV 


(To,e2//i 


581 


1 


15 


1 • 1012 


2 - 1018 


15 


2 


0,1 


945 


1 


15 


1,5-1012 


1 • 1018 


14 


2,5 


0,03 


946 


1 


15 


1,7-1012 


1,2 - IOI8 


13 


1,5 


0,1 


484 


5 


15 


1,3-1012 


1 ■ 1018 


26 


3 


0,02 


485 


5 


15 


1,5-1012 


1,3 • IOI8 


13, 26 


2,5 


0,03 


200N1 


20 


15 


2 - 1012 


1,2 • IOI8 


25 


2 




200N2 


20 


15 


3 - 1012 


2 • 1018 


16 


InT 





The conductivity of heavily doped sample (200N2) at 
low temperatures weakly decreases with a temperature 
decrease which corresponds to weak localization regime 
and has been reported in'^.. Note that at high temper- 
atures the samples demonstrated activated behavior of 



The maximum observed at high temperatures is at- 
tributed to a competition between scattering by optical 
phonons (/i oc T~^) and scattering by ionized impurities 
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FIG. 2: Temperature dependences of the Hall mobility for 
several samples with £4 conduction. 

(/i oc T . At temperatures 40-10 K one observes a further 
increase of /i with a temperature decrease. Such a be- 
havior is not expected for the nearest neighbor hopping 
(£3) since in the latter case the corresponding "mobihty" 
exponentially decreases with a temperature decrease^. 

Now let us consider temperature behavior of the con- 
ductivity and Hall effect for these samples at high tem- 
peratures 20-300. At T = 300—50 there exist an activated 
behavior resulting from an activation of the holes from 
occupied impurity band to the valence band. The value of 
El it is natural to ascribe to a distance from ep to the per- 
colation level within the valence band. Since the degree of 
compensation is small, the width of the impurity band is 
also small and, as it will be shown below ej? — << £o- 
Thus at small enough temperatures £1 ~ £9. However 
with the temperature increase when the value of the holes 
concentration p starts to be of the order of the dopant 
concentration one expects a transition from the regime 
£1 ^ £0 to the regime £1 ^ £o/2 — 14 meV. It is this be- 
havior which is observed experimentally (Table 1). This 
fact evidences that the low temperature behavior can not 
be attributed to £2 conductivity. 

An important feature is related to the sign of the Hall 
coefficient. As it is clearly demonstrated (Fig. 3 ) the 
sign at 300 K is opposite to the sign at small tempera- 
tures. 
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FIG. 3: Hall voltage as a function of magnetic field for T = 
4:.2K and T = 300K for a fixed direction of the current 1 = 2 
nA 

At the same time the concentration estimated from the 



Hall coefficient appears to be 2-3 orders of magnitude 
lower than at T = 300K. 

A specific attention was paid to the nonlinear effects 
in conductivity. Fig. 4 shows temperature behavior of 
resistance for different samples measured in a regime of 
constant current (0,1-2 nA). 
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FIG. 4: Temperature dependences of resistance for different 
samples measured in a regime of constant current (0,1 - 1 nA). 

It is seen that at low temperatures the samples ex- 
hibit a sharp transition from insulating behavior (Arre- 
nius law) to a metallic one; a magnitude of the resis- 
tance is decreased by 1.5-2 orders of magnitude. We 
believe that transition results from electric breakdown 
which takes place when the voltage drop on the sample 
(increasing with a temperature decrease) exceeds some 
critical value. The temperature of the transition differed 
for different samples and was equal to 3, 2,4, 2,1. This 
parameter is supposedly controlled by the voltage drop 
at the given temperature which depends on the slop of 
temperature behavior of the conductivity. 

Fig. 5 shows TV curves of the samples (obtained in 
the regime of constant /) at different temperatures. The 
curves clearly demonstrate a presence of S-shaped re- 
gions. One notes that the breakdown takes place at rela- 
tively weak electric fields which can be estimated as less 
than 10 V/cm. 
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FIG. 5: I-V curve for sample 485 at 4,2 K. 

This breakdown can hardly take place in the regime of 
hopping conductivity. The hopping conductivity in finite 
electric fields was studied earlier for the similar group of 
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samples^- . As it has been shown, the hopping conduc- 
tivity strongly (exponentially) increases with an electric 
filed increase, however this increase is still a gradual one 
and no breakdown-like behavior is observed (which agrees 
with theoretical predictions - see e.g.—) 

The breakdown can naturally result from an impact 
ionization of the acceptors to the valence band. However 
in this case the conductivity is expected to increase much 
larger than is observed experimentally. Then, for our case 
of relatively deep localized states, the breakdown would 
be expected at the fields at least orders of magnitude 
stronger than is observed experimentally (see e.gi^. 



III. DISCUSSION 



As it was noted earlier, the observed Hall mobility in- 
creases with decreasing temperature which excludes £3 
hopping conductivity as a mechanism leading to acti- 
vated behavior of conductance at low temperatures and 
leads us to a conclusion that the activated minority car- 
riers are delocalized. In this case the increase in the mo- 
bility with decreasing temperature can be explained by 
the scattering of delocalized carriers by acoustic phonons. 
This assumption does not contradict the fact that the 
mobility increases with temperature at high tempera- 
tures. We attribute this increase to the contribution of 
carriers activated to the valence band, where scattering 
on charged impurities prevails at corresponding temper- 
atures. As known, scattering on acoustic phonons is very 
sensitive to the effective mass of the carriers. On the 
other hand, the effective mass of delocalized carriers in 
the impurity band can noticeably exceed the carrier mass 
in the valence band as, in particular, was shown in our 
paper [3], where this excess was estimated by a factor 
of 2-3. For this reason, the changes in the scattering 
mechanism and a character of temperature dependence 
of mobility can be expected with decreasing temperature. 

There is another important evidence to support our 
picture which is related to the fact that a sign of Hall 
effect, observed at 4.2 K, is opposite to the sign of the 
Hall effect at 300K. The change of the sign apparently 
takes place at temperatures 20 - 40 K while the Hall 
data at these temperatures are not completely reliable. 

Indeed, with an assumption of an existence of a band 
of delocalized states(BDS) within the narrow impurity 
band we conclude that there are 2 mobility edges sep- 
arating the delocalized states from localized tail states. 
At low temperatures the conducting electrons are mainly 
concentrated near the bottom of the BDS where the ef- 
fective mass of the electrons is expected to be positive 
since the corresponding density of states decreases with 
energy decrease. 

Thus the Hall effect at high temperatures is produced 
by holes with positive effective mass while the Hall effect 
at low temperatures - by electrons with positive effective 
mass which explains the observed behavior. 



Note that the effective mass of an electron near the 
opposite mobility edge (i.e. near the top of BDS) is ex- 
pected to be negative. One has in mind that at high 
enough temperatures 20-50 K (which are larger than the 
activation energy £4) all the electrons are expected to be 
distributed over the BDS with nearly equal probability 
and thus the net Hall voltage is expected to be small due 
to different signs of the effective mass at different parts 
of the BDS. It gives an alternative explanation of the be- 
havior of Hall mobility which nearly vanishes at 20 -40 K 
but increases with temperature decrease when the con- 
tribution of electrons with positive effective mass start to 
dominate. 

An important evidence supporting our scenario is given 
by the character of non-Ohmic behavior. The fact that 
the conductance in low-Ohmic state is still much less than 
high temperature conductance over the valence band ex- 
cludes a breakdown to the valence band. Then, the elec- 
tric fields in our experiments are orders of magnitude less 
than necessary for impact ionization of deep centers^. 

In addition, extremely small electric powers (note that 
the breakdown behavior is observed at currents less than 
1 nA!) do not allow any heating effects. 

Thus we believe that the behavior observed results 
from the impact ionization of the electrons localized in 
the tail of the impurity band to the region of delocalized 
states. The important factor is a presence of a gap be- 
tween the Fermi level at and the "mobility edge" within 
the impurity band separating strongly and weakly local- 
ized states. This gap is filled by "intermediate" localized 
states (see fig. 6) with a concentration much larger than 
the electron concentration. 
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FIG. 6: Scheme of the density of states in the samples with 
narrow impurity band, where ei and £4 are the activation 
energies, eo is the isolated acceptor position, e/ is the Fermi 
energy, and Ad is the width of the band of delocalized states. 



The electrons existing within the BDS (responsible for 
the impact ionization) with much larger probability can 
be trapped by these "intermediate" states than by their 
initial centers. Since electron transport between the lo- 
calized states is a hopping one, the return of the electrons 
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to their initial position within the tail is suppressed. In 
its turn, the energy positions of these intermediate states 
are closer to the mobility edge than the tail states. As 
the result, the threshold field for the impact ionization 
from these intermediate states is lower that the threshold 
field for the ionization of the tails states. Correspond- 
ingly, the low-Ohmic state characterized by a presence 
of a finite electron concentration within the BDS can be 
supported at the electric fields which are lower than nec- 
essary for the ionization from the deep tail states which 
explains instabilities and S-shaped IV curves. (Note that 
the mechanism of a formation of S-shaped IV-curves re- 
lated to a presence of "intermediate" states between the 
ionized level and the conductance band was considered 

inia) 

The detailed theoretical treatment supporting the sce- 
nario given above is given in the Appendix. We have 
taken into account both kinetics of impact ionization and 
nonequilibrium electron distribution in strong enough 
fields. The first critical field, according to our calcula- 
tions, is given by the estimate 

ei?i(/Je-p.))^/'=^^ (1) 

Here U = VTi is an electron mean free path (r^ is transport 
relaxation time) while le-ph — vr^-ph where r^^ph is an 
electron-phonon relaxation time (note that both U and 
le-ph correspond to electron energies equal to £4). This 
equation has a clear physical explanation. The length 
{lile-phY^"^ is a length of electron diffusion during its in- 
elastic relaxation time (controlled by phonons) . Thus the 
critical field Ei allows the energy gain from the electric 
field of the order of the energy £4, the latter is just the 
energy necessary to activate an electron from the Fermi 
level to the delocalized states. 

The critical field E2 is given by the similar equation, 
however the energy £4 is in this case replaced by smaller 
energy £g corresponding to a typical distance from the 
electron trap states to the mobility edge. Thus we clearly 
have El > E2. 

Now let us give a rough estimates of the parameters 
allowing the scenario mentioned above. Our sample ex- 
hibits an increase of mobility with temperature decrease 
nearly saturated at T ~ 10 — 15 K at the values of 
the order of 10^T^~^cto^s~^. The saturation demon- 
strates a change of scattering mechanism (supposedly 
from electron-phonon scattering to electron-defect scat- 
tering). Implying a value of m ^ 10"^^g (which is sev- 
eral times larger than for the valence band) one estimates 
the transport relaxation time as 10~^^ s. To estimate 
Te-ph,s one concludes that transport relaxation rate 
at the crossover temperature only by a factor 1 /2 is con- 
tributed by phonons. If we assume that characteristic 
electron energy is some less than this saturation temper- 
ature like ^ 1 meV, one concludes that Te-ph,s ^ 4-10^^^ 
s since rjl^^^ ^ cx e^. Then, for the energy ~ ~ ImeV 
one estimates the electron velocity as ~ 3T0®cm/s . Thus 
the critical field Ei can be estimated as of the order of 



IQV/ cm. The electric field E2 is expected to be several 
times smaller, first, due to the fact that £g < e^, then 
also due to a possible role of electron-electron scattering 
which can be effective at large enough concentration of 
the mobile carriers n (see Appendix). 

As it is seen from Fig. 5 , our rough estimates of the 
threshold field are in order of magnitude agreement with 
experimental results. 

Then, from Fig. 3 it is seen that the conductance of 
the low-Ohmic state is of the order of the Ohmic con- 
ductance at r ^ 10 K and is somewhat smaller than 
the conductance at intermediate temperatures 20 - 40 K 
(when the electrons are effectively activated to the BDS). 
This fact evidences, first, that we do not deal with impact 
ionization of the holes to the valence band (the resulting 
conductance would be of the order of the high temper- 
ature one which is several order of magnitude higher). 
Then, it shows that only a fraction of localized electrons 
is ionized to the BDS. 

Note that the estimated threshold field is surprisingly 
low. We know only few references reporting S-shaped 
I-V curves at such a low fields, but for the case of shal- 
low impurities [ ]. Since in our case we deal with rel- 
atively deep acceptors we can conclude that the obser- 
vation of S-shaped I-V curves at low fields evidences a 
vicinity of metal-insulator transition which can not be 
of the Mott type. Indeed, it would contradict to the 
high-temperature activation energies which almost coin- 
cide with ionization energies of isolated acceptor while in 
a situation of the Mott transition Fermi level would be 
situated between energies of single occupied and doubly 
occupied states. In combination with the previous data 
it gives a strong support to the scenario of virtual An- 
derson transition. Then, it also allows to conclude that 
the observed £4 conductivity can not be attributed to £3 
channel. 

Thus, the results indicate that activation occurs to the 
band of the delocalized states appearing due to the An- 
derson transition in the impurity band. It is worth noting 
that the Anderson transition is typically considered in the 
single-particle picture in the absence of electron-electron 
correlations. In our case, the repulsion of holes from ac- 
ceptors prevents transport of delocalized holes even if the 
wave functions in single-particle approximation are delo- 
calized. However, owing to the activation of electrons 
from the Fermi level located at the tail of the impurity 
band to the indicated states, transport through such de- 
localized states becomes possible. In contrast to the real 
metal-insulator transition (when the Fermi level is lo- 
cated in the delocalized-state region), the manifestation 
of the delocalized states in our case can be described 
as the virtual Anderson transition. As we have shown 
earlier—, in the case of narrow impurity band existing in 
non-compensated samples due to the weak disorder po- 
tential the critical concentration of such a virtual Ander- 
son transition is smaller than given by the Mott criterion 
by the factor 

ln-2(eo/Ae) (2) 
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where Ae << Eq is a width of the impurity band. Cor- 
respondingly, although the dopant concentration in our 
samples was less than the critical one for the real Mott- 
Anderson transition, it can readily be larger than the 
critical concentration for the virtual Anderson transition. 

There is a peculiar question concerning localization ac- 
cording to Lifshits scenario^^. Indeed, the scatter in inter 
impurities distances can impose localization as well as a 
scatter of energies and the corresponding criterion is ex- 
cept a numerical coefficient is similar to Mott criterion. 
However, first, the coefficient for the Lifshits transition 
is still unknown and we can expect that the Lifshits lo- 
calization starts to be ineffective at lower concentrations 
than Anderson localization. Then, any deviation from 
purely random distribution of impurities (which can oc- 
cur in course of doping procedure) suppress Lifshits lo- 
calization. 

To conclude, we have given experimental and the- 
oretical evidences of an existence of virtual Anderson 
transition, characterized by an appearance of delocalized 
states within an impurity band above the chemical po- 
tential. Such a transition is expected to be typical for 
strongly doped, but uncompensated materials and pre- 
cedes (as a function of dopant concentration) standard 
Mott-Anderson transition. In our experimental studies 
of highly doped uncompensated p-type layers within the 
central part of GaAs/AlGaAs quantum wells it mani- 
fested itself, in particular, in activated behavior of low 
temperature conductivity. We have shown that the lat- 
ter can not be associated with known S2 or £3 mechanisms 
but results from an activation of electrons from the chem- 
ical potential to the delocalized states mentioned above. 
The fact that the type of the carriers in this case is dif- 
ferent from the one supporting conductivity within the 
valence band is evidenced by a change of the sign of Hall 
coefficient with a temperature decrease. Another evi- 
dence of the virtual Anderson transition is given by elec- 
tric breakdown observed at low temperature in unusually 
weak (for relatively deep centers) electric fields. We be- 
lieve that this behavior is related to impact ionization of 
the localized electrons from below the chemical potential 
to the band of Anderson-delocalized states. 
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V. APPENDIX 

We consider impact ionization of localized " minority" 
carriers - electrons - situated deep in the tail of the 
impurity band (with a concentration A'^i, controlled by 
compensating defects) by delocalized minority carriers 



to delocalized states with energies higher than the mo- 
bility edge Em (separating strongly and weakly localized 
states). We also have into account that in between 
and the position of the chemical potential fi at the equi- 
librium there exists a wide region of localized states which 
at low temperatures are occupied by holes and thus they 
can trap mobile electrons; the total number of such local- 
ized states N2 is much larger than iVi . While the analyt- 
ical solution of the problem with continuous spectrum of 
the localized states seems to be impossible, we will model 
our system as the 3-level one, including the "tail" states, 
the intermediate level of localized states separated from 
delocalized states by some energy Sg < Em, and the band 
of delocalized states. 

Denoting the the number of electrons captured on the 
localized states as n we describe the process of impact 
ionization at small temperatures as follows: 

dn , . . . , 

— = -BRn{N2 + n) + AinT2n + AinTi{Ni ~ h - n) 

dn ~ 

— = —Bj^n{n + fi) — AjJ-2nh + Bj^n(N2 — ^)(3) 
at 

{Bji is the coefficient describing a recombination of the 
electrons to any localized states, Bb. is the coefficient 
describing the recombination of electrons from interme- 
diate states to their initial positions, Aj is a coefficient of 
impact ionization, Ti and J-2 are the coefficients describ- 
ing the relative numbers of mobile electrons with kinetic 
energies larger than Sm and Eg, respectively. 

Note that in principle the coefficients Aj, Bji are en- 
ergy dependent. However we will neglect such a depen- 
dence since the width of the impurity band is not much 
larger than Em and thus the relative change of e at en- 
ergies higher higher than the threshold value £,„ is not 
large. 

Analyzing the first of EqsO one concludes that since 
N2 >> Ni if initially n = 0, n = 0, this solution is stable 
with respect to small fluctuations up to threshold electric 
field El corresponding to 

AiTiNi = BrN2 (4) 

At the same time for flnite n, n the solution of Eq|3]gives 

- BrN2 

Ni-fi (5) 

Here we make use of the relations A1T2 » Br (which 
holds at electric flelds of the order oi Ei Br« Br, 
N2 » Ni. This solution is stable until ii < Ni which 
holds until 

A1T2N1 > BrN2 (6) 

which defines another critical field £2- Since Em > Eg it is 
expected that for a given field T2 » Ti. Thus the value 
of E2 can me smaller than Ei and thus we deal with two 
branches of the solution: the one corresponding to n = 0, 
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n = which is stable with respect to small fluctuations at 
E < El and the solution given by EqO For electric fields 
E^Ei'ii corresponds to practically complete ionization 
of the electrons to the delocalized states. However this 
solution can exist at smaller fields until they are larger 
than E2, that is E2 < E < Ei the two solutions coexist 
and are stable with respect to small fluctuations. We will 
give more detailed analysis later with an account of the 
fact that both the function T and coefficients Bji can be 
different for the two branches. 

Now let us calculate the principal functions J-i and 
J-2- Considering the delocalized electrons we will take 
the mobility edge as the origin for e. By definition 

^1 = n-' f dei^ie)f{e), = n'^ f dsiy{e)f{e) 

(7) 

where / as an electron distribution function. 

We would fike to emphasize that, strictly speaking, 
the electron transport in our situation can hardly be de- 
scribed in a same way as for standard energy band - at 
least in the vicinity of the mobility edge, so our simple 
equations can be considered as semiqualitative. For rel- 
atively small n (when electron-electron processes are not 
effective) we will use a concept of energy diffusion and 
write with a neglect of electron-electron scattering 



Te-ph.T 



"^e—ph^s 



(8) 



Here the first term describes a diffusion of an elec- 
tron along energy axis due to interactions with thermal 
phonons with energies hwT = T (phonon absorption and 
stimulated emission characterized by the relaxation time 
Te-ph,T) and due to energy gains and energy losses in 
course of chaotic motion in electric field. 



(eEvTi 



(9) 



Ti is momentum relaxation time. The second term de- 
scribes an energy drift along energy axis due to sponta- 
neous emissions of phonons with typical frequency hujs 
characterized by relaxation time Te-ph,s- The electron- 
phonon relaxation for 2D electrons it is affected by the 
fact that the normal component of the phonon momen- 
tum is not conserved^^. According tc^^ at low temper- 
atures tJ1\j^ oc T^^e^^/^ which is similar to electron- 
phonon relaxation time in 3D metals and Hujt ^ T. 
At the same time for relatively small energies e < 
{ms'^Wy^^ (where s is a sound velocity while W is an en- 
ergy of lateral quantization which in our case of 2D impu- 



rity band is of the order of the Bohr energy) rjjp^ s ^ 
while hus ^ e. 

The solution of Eql8] depends on the relation between 
De and phonon contribution to the energy diffusion. If 
the last term dominates, the solution corresponds to equi- 
librium, / oc exp— (e/T). If we have the opposite rela- 



tion, which holds when 

we deal with strongly nonequilibrium distribution 

e 



f cx exp — / de' 
^0 



[eEYv'^TiTe-ph,s 



(10) 



(11) 



According to our experimental data on temperature de- 
pendence of the mobility, at the energy region of interest 
e ^ £m Ti oc e^"^ where d ~ 2 and thus 



/ oc exp - 



5ieE)^v^nie)Tsie) 



This function is nearly constant at e < ei where 



5{eEfv'^T,Te-ph^. 



= 1 



(12) 



(13) 



and extremely steeply (much stronger than simple expo- 
nential) decays at e > ei. As it is seen, for strong enough 
fields when £1 >> T the inequality of EqlTU] holds auto- 
matically. 

For the further analysis we should have estimates for 
Bn and Aj. To estimate Aj we have in mind that Ai = 
vaj where aj is the ionization cross-section. The lat- 
ter (in 2D) can be very roughly estimated as ac{(Jt/^'^) 
where ac ~ {e^ /{Ke)y^^ is the geometrical cross-section 
of a Coulomb scattering of an electron with kinetic energy 
e > Em by a charged center, k being a dielectric constant, 
at is a cross-section of an elastic scattering of an electron 
by the trap center and A is a typical electron wavelength 
which can be roughly estimated as ~ N3 where 
is a concentration of the delocalized states. Note that 
at can be estimated from above basing on the known 
values of electron mean free path k as at ~ (iV2^i)^^ 
(actually the corresponding cross-section can be smaller 
if electron momentum relaxation is dominated by other 
mechanism). As for we shall consider recombination 
due to phonon emission. In this case recombination pro- 
cess involves simultaneous interaction of an electron with 
phonon and trap center (neutral). One can estimate the 
recombination cross-section an Bji/v as at /vTe-ph,s- 

Correspondingly, the estimates of the critical fields Ei 
{E2) can be rewritten as 



BrN2 



AiTi 



No 



(2) 



N2,NiVTe-ph,s<yC Pl{2) 



1 



(14) 



Assuming that ~ -^3, one concludes that our scenario 
can hold provided 



(15) 



since in any case J- < 1. As it can be readily estimated, 
ac ^ 2 • 10~^cm while VTf,-ph,s ~ 10~^ cm. Thus this 
criterion is obeyed for A^i > 0.5 ■ 10^ cm~2. 



8 



Since an equilibrium phase at T — > corresponds to 
n = 0, the threshold field and, correspondingly, is 
related to a distribution given by Eq[T2] which can be 
modelled as a cut-off of of /(e) at e = ei. Correspond- 
ingly, Tx ^ - em)(ei - £,„)/ei. Thus Ex can be 
estimated from an equality £i ~ which gives 



1/2 



5172 



(16) 



where = ivTi)]^^^^, le-ph,m = {vTe-ph,s)\e=6,^- In- 
deed, for E > El 1 and the condition of EgfTH 
holds. Note that by definition coincides with the ac- 
tivation energy £4. 

As it is seen, the estimate of E2 is a similar one and 
we have E2 < Ex because of the replacement of £„ by 



£f7 ^ £r) 



Eqdni 



If n is high enough, the energy relaxation is dominated 
by electron-electron rather than by electron-phonon pro- 
cesses and thus 



/ = exp -e/Te 



(17) 



However the energy transfer from the electron system to 
the bath is still due to electron-phonon processes. The 
effective electron temperature Tg which is established in 
the system can be derived from the energy balance be- 
tween energy gain from electric field aE^ = n(eE)'^Ti/ra 
and its decay to the thermal bath ~ n(Tf, — T)/Te-ph. 
Thus we have 



eHt.) = ^ - 1 



(18) 



As it is seen from EqlT7]and EgfTSl formally the estimate 
for Te except of a numerical factor coincides with esti- 
mate for £1 and characterizes an average electron energy 
which is controlled by energy balance between electrons 
and the lattice and does not depend on the efficiency 
of electron-electron processes. However the asymptotic 
behavior at large energies (> Te,£i ) strongly depend 
on these processes. Namely, without these processes the 
decay of f{e) is significantly stronger within the region 

£1 < £. 

Let us assume that for electron concentrations n = Nx 
electron-electron processes are much more effective for 



energy relaxation than electron-phonon ones. Thus the 
second branch (where n ^ Nx) corresponds to effec- 
tive electron-electron scattering and electron distribu- 
tion is described by Eq[T71 Correspondingly, T2 = 
exp(— £g/re). Thus the electric field E2 is given by the 
equation 



Te{E2) = Sgln-^ 



AiNx, 

BrN2' 



(19) 



However for the second branch where electron-electron 
scattering is assumed to be effective, we shall also con- 
sider Auger processes. Following the same rough proce- 
dure as applied above, we can estimate the cross-section 
(jR as ~ {(JtYacn. Note that in this case Eql5] has a 
solution 



Nx 



l+iV2/(^3^2) 



(20) 



and thus Eq. [6] holds automatically. However one has in 
mind that the Auger processes dominate only if 



(21) 



which holds only for large enough n and is violated at 
some critical n = ric- At n < ric electron-phonon scatter- 
ing dominates and the concept of electron temperature 
is not valid, so the system can not stay at the second 
branch. Thus the estimate of Eq[T9]can be rewritten as 



Te{E2) 



1 -uN^Nx. 
£oln ( 



(22) 



Note that the electron-electron processes can lead to ad- 
ditional source of multistability. Indeed, the concentra- 
tion of high energy electrons given by Eq |17l can be at 
the same electric field larger than given by EqlT^ Cor- 
respondingly, the effective impact ionization by electron 
distribution of Eq. [17] can be supported at electric field 
somewhat weaker than necessary to support the distri- 
bution of EqlH 
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